Phase diagram of a Heisenberg spin-Peierls model with quantum phonons 
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Using a new version of the density-matrix renormalization group we determine the phase diagram 
of a model of an antiferromagnetic Heisenberg spin chain where the spins interact with quantum 
phonons. A quantum phase transition from a gapless spin-fluid state to a gapped dimerized phase 
occurs at a non-zero value of the spin-phonon coupling. The transition is in the same universality 
class as that of a frustrated spin chain, to which the model maps in the anti-adiabatic limit. We 
argue that realistic modeling of known spin-Peierls materials should include the effects of quantum 
phonons. 



PACS numbers: 75.10.Jm, 75.50.Ee, 71.38. +i, 71.45.Lr, 63.20.Kr 

Challenged by the discovery of high-temperature su- 
perconductivity in doped antiferromagnets, our under- 
standing of quantum magnetism in low dimensions has 
increased significantly over the past decade jy. However, 
the effect of the interaction of quantum spin systems with 
further degrees of freedom such as disorder, phonons, and 
holes produced by doping is still poorly understood. In- 
terest in models of spins interacting with phonons has in- 
creased significantly since the discovery of a spin-Peierls 
transition in the inorganic compound CuGe03 [Q. The 
availability of large, high-quality single crystals has led 
to much more extensive experimental studies |3| than on 
the organic spin-Peierls materials studied in the 1970's 

The fact that a spin- 1/2 antiferromagnetic Heisenberg 
chain is unstable to a static uniform dimerization is 
known as the spin-Peierls instability. This occurs because 
dimerization opens a gap A in the spin excitation spec- 
trum and lowers the total magnetic energy at a greater 
rate than the increase in elastic energy due to the dimer- 
ization. Until very recently almost all theoretical treat- 
ments have used this static picture which we can only 
expect to be valid in the adiabatic regime where typical 
frequencies u> of the phonons associated with the dimer- 
ization are much smaller than the magnetic energy scales 
such as A and the antiferromagnetic exchange integral J. 
It has recently been pointed out that CuGeC>3 is not in 
this adiabatic regime |^^,|| , stimulating several numeri- 
cal studies with dynamical phonons [^|^] . 

In this Letter we study a model of a spin-1/2 antifer- 
romagnetic Heisenberg chain interacting with quantum 
phonons using a powerful new numerical technique that 
allows an essentially exact treatment of both the spins and 
the phonons at a fully quantum-mechanical level. Our 
main result is the phase diagram in Fig. [I] in which the 
adiabaticity parameter J/ui varies over several decades. 
We find that the spin-phonon coupling must be larger 
than some non-zero critical value for the spin-Peierls in- 
stability to occur. This is in contrast to the static case 



(u>/ J — > 0) for which dimerization occurs for any value 
of the coupling. Hence, quantum lattice fluctuations can 
destroy Heisenberg spin-Peierls order. We find that the 
quantum phase transition from the spin-fluid state to 
the gapped state is in the same universality class as the 
dimerization transition of the J\-Ji frustrated spin chain. 
Our results have important implications for the modeling 
of spin-Peierls materials. 

The model we study is one of the simplest possible. It 
consists of a local phonon on each site and the antiferro- 
magnetic exchange on neighbouring sites varies linearly 
with the difference between the phonon amplitudes on 
the two sites. The Hamiltonian is 
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H = J2[J + 9 [bi+i + 6.1+1 - bi - 4)) §i ■ S l+1 + u ]T b\b 



(1) 



Here Si is the 5 = 1/2 spin operator on site i and hi 
destroys a phonon of frequency u> on site i. We assume a 
periodic chain of N sites. 

Insight into this model can be obtained by consider- 
ing the anti-adiabatic limit (w >> J). One can then 
integrate out the phonon degrees of freedom to obtain 
the following effective Hamiltonian for the spin degrees 
of freedom |l(| 



N 



N 



7~L e fi = <h } Si ■ Si. 



J2 J]] Si ■ S; 



i+2 



(2) 



where J± = J + g 2 /uj and J2 = g 2 /2u>. Uhrig [|| recently 
obtained the same Hamiltonian, calculating J\ and 3i to 
next order in J/lu. 
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J + g 2 /uj - 3g 2 J/2uj 2 
g 2 /2uj + 3 9 2 J/2lu 2 + . 
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The frustrated spin chain Eqn. (|^) or J1-J2 model has 
been extensively studied and is well understood. If a = 



1 



Ji] J\ then at a critical value of a = a c — 0.241167(5) 
the model undergoes a quantum phase transition from a 
gapless spin-fluid state with quasi-long-range antiferro- 
magnetic order to a gapped phase with long-range dimer 
order 11 Uhrig pointed out that this implies that in 
the anti-adiabatic regime (Q) possesses a non-zero critical 
coupling g c . To second order in J/lu, 



a c J 



l/2-a c + 3(l + a c )J/2w' 



(5) 



We have confirmed this result numerically (see Fig. [I]). 
Furthermore, this non-zero critical coupling g c still oc- 
curs well into the adiabatic regime. It is interesting that 
although (H) is only valid to second order in J/lu it gives 
a good description of g c up to J/ui ~ 1 . 

Models such as (|l|), which involve bosons are a chal- 
lenge to study numerically due to the large number of 
degrees of freedom per site. The density matrix renor- 
malization group (DMRG) method |l3| has the poten- 
tial for obtaining definitive results for these models by 
studying very large systems. Several schemes based on 
the DMRG have recently been developed to treat models 
involving phonons p4] , p^|Jl6| . We employ a new "four- 
block" DMRG method ]16[| which allows us to treat the 
phonons and spins on an equal footing and to study sys- 
tems as large as 256 sites. This is in contrast to some re- 
cent exact diagonalization studies of spin-phonon models 
that were limited to small systems and/or used uncon- 
trolled truncations of the phonon degrees of freedom §,§. 
We previously used this method to obtain the phase di- 



agram of the Holstein model with spinless fermions |17|. 

The four-block method can be used to calculate the 
ground state energy E and the singlet and triplet gaps 
A ss and A st for periodic systems [jl6[ . Tabic | shows the 
DMRG convergence of the gaps with the single trunca- 
tion parameter e 1 18 for a representative parameter set. 
It can be seen that the gaps are sufficiently well resolved 
to be useful for finite-size scaling analyses. The error of 
around 0.1% in the A = 128 site system is typical of 
the error in the largest systems studied for a given set of 
parameters. 

We determine the critical coupling using the gap- 
crossing method used by Okamoto and Nomura ]l^ ] 
to determine the critical coupling a c in the frustrated 
Heisenberg model (|2j). The convergence of the crossover 
coupling a c (N) with A is rapid due to the absence of 
logarithmic corrections at the critical point |Tl],[l2|,[l9| . If 
the system is gapless with quasi-long-range Neel order 
for < g < g Cl the lowest excitation is the triplet state, 
i.e. A st < A ss (for sufficiently large A) and A st , A ss — > 
as A — > oo. If for g > g c the system has a non-zero 
gap A and is dimerised with a doubly degenerate ground 
state, then the first excited singlet state becomes degen- 
erate with the ground state in the bulk limit That 
is, A ss < A st (for sufficiently large N), A ss — * 0, and 



A s t - > A > as A ^ oo. A finite lattice crossover 
coupling g c (N) is defined by A st = A ss . As shown in 
Table |l[ g c (N) rapidly approaches a limit as A — > oo. 
This limit is the critical coupling g c separating gapless 
and gapped phases. For the J/lu > 1 cases, where the 
A dependence is substantial, g c {N) is well described by 
the functional form g c (N) ~ g c — Aexp(-BN) and non- 
linear fitting is used to determine g c pOfl . The resulting 
phase boundary is plotted in Fig. |l| The DMRG, dis- 
cretization and fitting errors in g c are estimated to be no 
greater than a few percent. 

From conformal invariance the finite-size energies of 
the spin-fluid should satisfy E3]: 
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where is the bulk ground state energy density and 
Vo = Vi = v a is the spin wave velocity. The combination 
of the gaps in Eqn. (Q) is chosen to cancel the logarithmic 
corrections. 

We have performed a number of consistency checks on 
our results. First, in Fig. ||, v a as determined by (||) is 
plotted as a function of g for a phonon frequency deep 
into the anti-adiabatic regime (J/lu = 0.005), together 
with the same quantity determined from the correspond- 
ing J1-J2 model. The results from the two approaches 
agree well. This confirms the mapping between the two 
models in the anti-adiabatic regime. Second, we note 
that the DMRG results for the phase boundary agree 
well with the result (^J) from the mapping in the anti- 
adiabatic limit (See the dotted line in Fig. [l]). Third, for 
general phonon frequencies, we calculate the ratio vq/vi 
which should equal unity. At g = g c it is one within 
errors expected from corrections to scaling and DMRG 
truncation, over the range of frequencies studied. Values 
vary from 0.98 ± 0.04 for J/lu = 0.005 to 1.07 ± 0.10 for 
J/lu = 10. 

For a K-T transition, the gap A = liniAr^oo A st is ex- 
pected to have an essential singularity at g — g c . In Fig. 
|§|, A st is plotted as a function of g for various A in a 
case of intermediate coupling J/lu = 1 . Two-point linear 
extrapolations (in 1/A) to A = 00 are included in the 
plot. These estimates of A are shown to be well fitted 
by the K-T form p| A ~ Af(g) exp(-S(/( 5 )) 2 ) where 
f(g) = (g - 5c )- 1 / ? ~Note that the gap crossover method 
(Table is substantially more accurate than this fit- 
ting procedure for determining g Cl the latter tending to 
overestimate g c Em. 

In the adiabatic regime (lu << J) there is strong mix- 
ing between spin singlet and phonon excitations. An 
analogous effect was observed for the Holstein model [[l7| . 
In the case of (Q) this is manifest in nonlinear corrections 
to the scaling of A ss . That is, A ss is found to be phonon 
like (flat in 1/A) until the characteristic spin energy 
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2irJ/N decreases below the bare phonon frequency cu, at 
which point A ss begins to vanish, as X/N (0 < g < g c ), or 
exponentially (g > g c ). This can be seen in Table |] from 
the slow convergence of g c {N) with N for the J/u> = 10 
case. 

Next, we consider the validity of the static approxima- 
tion in the adiabatic regime, where the phonon operators 
bi in (1) are replaced by the constant dimerization (— l) l 5, 
the total energy is minimized as a function of S then the 
gap is calculated for this optimal value of S. This cal- 
culation was performed by using the four-block DMRG 
method to solve for the ground state energy and gap in 
the dimerised Heisenberg mode l |22|| . The resulting adia- 
batic curve is compared in Fig. |3j to the extrapolated gap 
A found from our numerical results for J/oj — 10. We see 
that even in this adiabatic region treating the phonons in 
the mean-field approximation is not fully reliable, partic- 
ularly for the purposes of quantitatively fitting the cou- 
pling g from the experimental triplet gap. The situtation 
is far worse for phonon frquencies relevant to CuGeC>3. 
For example, for the J/oj — 1 case in Fig. 0, the adiabatic 
curve would not fit on the same scale as the curve from 
the fully dynamical model. 

To consider our results in the context of experiment, 
estimates of a number of parameters for various spin- 
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Peierls compounds are listed in Table III. It can be seen 



from these estimates and our results that the static ap- 
proximation is highly questionable for CuGe03, and may 
not be valid for the organic spin-Peierls materials. A re- 
lated question is the use of an explicit next-neighbour 
(J2, frustration) term in adiabatic spin-phonon models 
ofCuGe0 3 " 



|21[ |. The value of J2 required to achieve 
agreement with susceptibility and magnetic specific heat 
data is generally very large (J2/J1 ~ 0.3). Attempts have 
been made to justify the inclusion of a J2 term on the 
basis of Cu-O-O-Cu superexchange paths ^1). However, 
Ref. Q and the present analysis suggests that an explicit 
J2 term may not be required in order to describe exper- 
imental results if the phonons are treated quantum me- 
chanically since the phonons induce a next-nearest neigh- 
bour interaction. 

To conclude, we have numerically determined the 
phase diagram of a spin-Peierls model (^) with high ac- 
curacy. Our results are consistent with a mapping of 
the model to the frustrated spin chain (|J) in the anti- 
adiabatic limit (large phonon frequency). For a wide 
range of phonon frequencies compared to the exchange 
there is a phase transition at a non-zero value of the cou- 
pling g from a gapless spin-fluid state to a gapped dimcr 
phase (2^| . The transition is in the same universality class 
as the Kosterlitz-Thouless transition in the frustrated an- 
tiferromagnetic chain (^) . Quantum phonon fluctuations 
are important in known spin-Peierls materials. 
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TABLE I. Four-block DMRG convergence of the singlet 
and triplet gaps A ss and A st of the spin-Peierls model (ji]) 
with the truncation parameter e for various periodic lattices 
of size N, where J/u = 1 and g/u = 0.4. 
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A ss /aj 


A s t/cj 
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10~ ib 


0.31374961 


0.5183251 
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0.31372889 


0.5183254 
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0.133925 


32 
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0.0765958 


0.133785 
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0.0765933 


0.133778 
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0.014909 


0.04009 


128 


lO" 11 


0.014817 


0.03856 


128 


IO" 13 


0.014619 


0.03790 


128 


1Q -14 


0.014648 


0.03775 



TABLE II. Convergence of the crossover coupling g c {N)/u) 
with lattice size N for various values of the adiabaticity pa- 
rameter J/u. g c {N) is defined by A ss = A s t and converges to 
the critical coupling g c as N — » oo. 
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0.3021 


0.274 
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0.223 


0.3087 
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0.3092 


0.318 


0.249 
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0.318 
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0.339 



TABLE III. Estimates of the exchange J, phonon fre- 
quency u, and energy gap A for various spin-Peierls materials. 
All are given in units of Kelvin. (We are unaware of any other 
measurements of the frequencies of the dimerization phonon 
in organic materials.) 
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FIG. 1. Zero temperature phase diagram of the spin-Peierls 
antiferromagnetic chain of spins interacting with quantum 
phonons (Eqn. ([j])). For small spin-phonon coupling g the 
system is a gapless spin-fluid. For large g the system is dimer- 
ized and has an energy gap. The diamonds with error bars 
denote the phase boundary from this DMRG study. The dot- 
ted line is (Eqn. (|E|)) the phase boundary which results from 
an approximate mapping onto the Ji- J2 model (frustrated an- 
tiferromagnetic chain, Eqn. (^)) which becomes exact in the 
anti-adiabatic limit J/u — > 0. 
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FIG. 2. The spin wave velocity v a as a function of the 
electron-phonon coupling g for a large phonon frequency 
J/u = 0.005 (anti-adiabatic regime). DMRG results (from 
extrapolating the N = 16 and 32 data using (^)) are indi- 
cated by diamonds. The triangles are obtained by solving the 
N — 16 and 32 site J1-J2 model (Eqn. (^)), with Ji and J2 
given by ([]) and (Q). The dotted line arises from applying an 
approximate result for v a pd] to the same J\-Ji model. 
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FIG. 3. The singlet-triplet gap A s t of the spin-Peierls 
model as a function of the coupling g for various lattice sizes 
N for an intermediate phonon frequency J/ui — 1. Extrapola- 
tions (in 1/N, using the two largest values of iV) to N = 00 are 
given by the solid diamonds. These are fitted to the K-T form 
Af(g)eM~Bf(g) 2 ), where f(g) = (g - g^ 1 ' 2 (solid line). 
The critical coupling g c is not obtained from this fit. It is 
substantially more accurate to use the gap crossover method 
(see Table nf). The inset shows the extrapolated gap (us- 
ing TV = 32 and 64) for a small phonon frequency (adiabatic 
regime) J/u = 10. The dashed line is the result for the static 
limit where the quantum phonon fluctuations are neglected 



